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ABSTRACT 

Polaris is one of the most observed stars in the night sky, with recorded observations spanning 
more than 200 years. From these observations, one can study the real-time evolution of Polaris via 
the secular rate of change of the pulsation period. However, the measurements of the rate of period 
change do not agree with predictions from state-of-the-art stellar evolution models. We show that this 
may imply that Polaris is currently losing mass at a rate of M ~ 10~ 6 Mq yr _1 based on the difference 
between modeled and observed rates of period change, consistent with pulsation-enhanced Cepheid 
mass loss. A relation between the rate of period change and mass loss has important implications 
for understanding stellar evolution and pulsation, and provides insight into the current Cepheid mass 
discrepancy. 

Subject headings: stars: individual (Polaris) — stars: evolution — stars: variables: Cepheids — stars: 
mass-loss 



1. INTRODUCTION 

Classical Cepheids are radially pulsating variable stars 
and are ideal laboratories for understanding stellar as- 
trophysics and evolution. Radial pulsation constrains 
stellar structure along with Cepheid masses and lumi- 
nosities via the Period-Mean Density relation. However, 
observations of the rate of pulsation period change di- 
rectly constrain stellar evolution because this measures 
the change of the stellar mean density as a function of 
ti me. 

iTurner et al.l (|2006l ) compiled measurements of rates 
of period change for a sample of Galactic Cepheids and 
showed that they broadly agree with rates of period 
change as a function of pulsation period computed from 
stellar evolution models. The authors also found that 
observed rates of period change for Polaris and DX Gem 
are larger than those of most Cepheids and are consistent 
with being a Cepheid on its first crossing of the instabil- 
ity strip. The first crossing occurs as a star evolves across 
the Hcrtzsprung gap after main sequence evolution on a 
Kelvin-Hclmholtz timescale. 

That work illustrated the power of measured rates of 
period change to constrain stellar evolution models, how- 
ever, one can go a step further. By combining observa- 
tions of the rate of period change with measurements 
of mean stellar radius, mean luminosity and/or effective 
temperature, we can directly test the input physics used 
in stellar evolution models. In this work, we produce a 
grid of state-of-the-art stellar evolution models and com- 
pare stellar evolution tracks with the measured radius 
and effective temperature of the first-overtone Cepheid 
Polaris. For evolutionary tracks that cross the Cepheid 
instability strip, we compute theoretical rates of period 
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change for comparison. From this comparison, we con- 
strain the physics of stellar evolution models and other 
fundamental properties of Polaris. 

Polaris is an ideal candidate for this analysis. It is an 
astrometric binary (as part of a multiple star system) 
with a dynamically measured mass of M = 4.5±?;t M Q 
(jKamperl 119961: lEvans et"aLl 120081 ). has been observed 
using K-band interferometry yiel ding a measuremen t 
of the its mean angular diameter ( Merand et al.l [2006f> . 
and h as a precise, revised Hipp arcos parallax measure- 
ment (|van Leeuwen et al.l I2007T) . The effective temper- 
ature of Polaris has also been measured spectroscopi- 
cally ([Usenko et al.l 120051) . The combination of these 
observations provide strong constraints for stellar evo- 
lution models. There have also been a number of mea- 
surements of the rate of period change for Polaris based 
on the combination o f rece nt and archival observations. 
ISprecklev fc Stevens! (|2008l ) measured a rate of period 
change P = 4.90 ± 0.26 s yr -1 based on 50 years of 
observations, while ITurner et al.l (|2005l ) found a rate of 
« 4.4 s yr^ 1 based on 160 years of observations. In 
this work, we present a rate of period change based on 
data from lTurner et al.l (|2005l ) with new data added from 
measurements over the past decade. 

In the next section, we describe the measurements of 
the rate of period change in greater detail. In Section 3, 
we describe the stellar evolution code and models used in 
this study, while in Section 4, we compare the predicted 
rates of period change with the observed rate where we 
find disagreement. In Section 5, we discuss the impli- 
cations of our study toward understanding the structure 
and evolution of Cepheids as well as the problem of the 
Cepheid mass discrepancy. 

2. OBSERVED RATE OF PERIOD CHANGE 

The recent comprehensive study by ITurner et al.l 
( 2005) of photometric and radial velocity measurements 
(up to 2004) of Polaris yielded rates of increase in the 
pulsation period of Polaris of P = 4.45 ± 0.03 and 
4.28±0.73 s yr _1 . This result is in good accord with pre- 
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Figure 1. O-C data from 1844/5-2011, from both radial velocity 
and photometry measures, are shown for the Cepheid Polaris. For 
the purposes of fitting, the data were divided into two sets (pre- 
1963 and post-1965) because of the apparent discontinuity in the 
O-C curve around 1963-66 (Epoch 2400-2600). The weighted, least 
squares quadratic fits to the data sets are also shown, where the 
blue solid line is the fit to the earlier dataset and the red dashed line 
is the fit to the later. The crosses represent data from Turner et al.l 
(2005) and circles our new timing data. The early data show a 
period increase of 4.44 ± 0.19 s yr — 1 and the later data show a 
period increase of 4.50 ± 1.45 s yr — x . 

vious period studies of Polaris. We added five photomet- 
rically determined t imings (covering 200 2-2011) to those 
given in Table 1 of iTurner etail (|2005f l. The resulting 
O-C analysis plot is given in Figure [1] As noted previ- 
ously, there is an apparent discontinuity in the O-C plot 
occurring near epoch 2400-2600. Analyses of the tim- 
ings prior to (1844/5- November, 1962) and after (July, 
1965-2011) this small anomaly yielded common values of 
period increase P = AAA ± 0.19 and 4.50 ± 1.45 s yr _1 , 
respectively (the larger uncertainty for the later period 
change is due to the more limited time span covered). 
When all data are taken as one set, the calculated pe- 
riod increase is 3.93 ± 0.08 s yr _1 . For the purpose of 
this study, we adopt the value of P = 4.47 ± 1.46 s yr _1 , 
as the average of the two fits on the two separate s ets de - 
scribed above. In fact, as noted bv lTurner et al.l (|2005l ). 
the discontinuity in the O-C curve coincides with the 
possible beginning of an accelerated decrease in the light 
amplitude of Polaris, however this requires further study. 

3. STELLAR EVOLUTION MODELS 

We co mpute state-of-the-art stellar evolution models 
using the lYoon fc Langen (12 005) version of the code de- 
veloped bv lHeger et al.l (|2000f ) . Models are computed for 
the mass range of 3 - 6.8 M & in steps of 0.1 M , con- 
sistent with t he me asured mass of Polaris determined by 
lEvans et al] (|2008ft . We also vary the amount of con- 
vective core overshooting. Convective core overshooting 
is parameterized as the distance that convective eddies 
penetrate above the core boundary in a star over its 
evolutionary timescale . The distance is defined, in the 
lYoon fc Langerl ()2005D code, to be a c Hp, where Hp is 
the pressure scale height and a c is a free parameter. In 
this work, we compute models for a c — 0-0.4 in steps of 
0.1. We vary the amount of overshooting because over- 
shooting during main sequence stellar evolution acts to 
produce a more massive helium core and hence a more 
luminous Cepheid for the same stellar mass with no over- 



shooting. There is no enhan ced Cepheid mass lo ss as- 
sumed like that explored by iNeilson et al.l (120111) and 
the m odels are assumed to have the iGrevesse fc Sauval 
(1998) solar comp o sition . Mass loss is treated using 
the Ide Jager et aLl (I1988D recipe for cool stars and the 
Kudritzk i et al.l ( |1989|) prescription for hot stars. 

Evolutionary tracks for models without convective core 
overshoot are shown in Fig. [2] for a sample of stel- 
lar masses along with the measured effective tempera- 
ture and luminosity of Polaris. The effec tive tempera- 
ture w as spectroscopically determined by lUsenko et al.l 
(pOOll to be T eff = 6015 ± 170 K, while the radius is 
R = 43 .5 ± 0.8 R P) based o n angular diameter measure- 
ments (iMerand et al.1 12006D and t he revised Hipparcos 
parallax (|van Leeuwen et ail I2007D . hence logL/L = 
3.347±0.061. In Fig. [21 not all models are consistent with 
the observed properties, suggesting the observations are 
precise enough to constrain models. In the next section, 
we compute the rates of period change from the stellar 
evolution models. 

4. RATES OF PERIOD CHANGE 

ITurner et al] (|2006l ) derived a relation for computing 
the rate of period change from the predicted stellar prop- 
erties of evolution models from the Period-Mean Density 
relation and n oting that the pulsation constant Q oc P 1 / 8 
(|Fernielll967D . 



P 
P 



6L 
7L 



24 7rff 

7 Teff' 



(1) 



This relation assumes that mass loss is insignificant. We 
compute rates of period change P using this relation 
where the values of L and T c g along with the time deriva- 
tives are given by the stellar evolution models. 

When the stellar evolution models predict radii and 
effective temperatures that fall within the boundaries of 
the Cepheid instability strip, we compute the model rel- 
ative rate of period change, P/P using Eq. 1. This is 
done to avoid making any assumptions of the pulsation 
period from the theoretical stellar models. We show the 
predicted rates of period change as a function of stellar 
luminosity in Fig. [3] The observed relative rate of period 
change is scaled by the calculated fundamental period for 
Polaris, Pq — 5.61 days. The calculated fun damental pe- 
riod is the overtone period plus 0.15 dex (jTurner et al.1 
12005] ). 

We find that evolutionary models with masses 5.4 - 
6.4 Mq and convective core overshooting a c < 0.4 pre- 
dict effective temperatures and radii consistent with ob- 
servations, but the predicted mass is degenerate with the 
assumed amount of overshooting. This can be seen in 
Fig. [2] for models with a c — 0, and masses, M = 5.9 - 
6.4 M Q ; smaller masses are measured for a c > 0. The 
de generac y between mass and a c is discussed in detail 
byLKeller (2008|). Increasing a c in stellar evolution mod- 
els decreases the blue loop width, however, this does not 
change the rate of period change for a given stellar lu- 
minosity during blue loop evolution. The rate of period 
change is dominated by the most important timescale 
at each stage of stellar evolution. On the Hertzsprung 
gap, a star undergoes thermal relaxation, meaning the 
timescale and the rate of period change is defined by the 
Kelvin timescale tk, which is oc M . The evolutionary 
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Figure 2. Stellar evolution tracks computed using the Yoon & Langcr (2005) code for initial masses ranging from 3 Mq to 6.8 Mq, 
assuming no convective core overshooting, a c = (left) and the region of the HR diagram containing th e effective temper ature and 
luminosity of Polaris(right). The black dashed lines represent the boundaries of the Cepheid instability strip (Bono ct al. 2000) while the 
boxes denote the effective temperature and luminosity of Polaris. 

timescale of a star on the blue loop is defined by the 
helium nuclear burning timescale, th c , which itself de- 
pends on the helium core mass Mh c , i.e. the luminosity. 
Thc helium core mass is determined by the combination 
of stellar mass and the amount of assumed convective 
core overshoot. Thus, the rate of period change for an 
evolution model with a given luminosity, mass, and no 
overshooting is approximately the same as for a model 
with the same luminosity, but a different amount of over- 
shooting and corresponding mass. 

The rates of period change are plotted in Fig. Has a 
function of luminosity for the models that have a pre- 
dicted radius and effective temperature consistent with 
the observed values during the blue loop evolution. This 
plot differs from the results in Fig. [3] such that the the- 
oretical value of the pulsation period is assumed to be 
that of Polaris. Any uncertainty introduced by this as- 
sumption will not affect the results. The stellar evolu- 
tion models predict rates of period change that occur in 
three regions for the plot in Fig. [3] The first region in- 
cludes rates of period change where P > 15 s yr -1 in 
Fig 2J which corresponds to stellar models rapidly evolv- 
ing across the Hertzsprung gap. The second region is 
for P < s yr _1 , that corresponds to blueward evolu- 
tion during the blue loop phase of stellar evolution, while 
the third region is < P < 1 s yr _1 corresponding to 
redward evolution during the blue loop phase. 

Thc observed rate of period change is inconsistent with 
any of the three phases of stellar evolution, and con- 
trary to p revious claims Polaris can not be a first-crossing 
Cepheid (iTurner et alj|2005l 12006ft This difference was 
also found bv lSprecklev fc Stevens! (|2008l ). The question 
is what physics is missing that can resolve this difference. 
We rewrite the rate of period change from Eq. Q] to be 



Figure 3. Predicted rates of period change scaled by pulsation 
period as a function of stellar luminosity computed from the stel- 
lar evolution models when the evolution tracks intersect assumed 
boundaries of the Cepheid instability strip (sec Bono ct al. 2000, 
for details). The different colors represent differing amount of as- 
sumed convective core overshooting, a c = (green), 0.1 (blue), 0.2 
(light blue), 0.3 (yellow), and 0.4 (magenta). The box represents 
the observed luminosity and the rate of period change range for 
Polaris. The inset zooms in on the luminosity range for Polaris 
and shows a distinct difference between theoretical and observed 
rates of period change. 
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Figure 4. Predicted rates of period change as a function of stel- 
lar luminosity computed from stellar evolution models when the 
models have effective temperature and radius consistent with ob- 
servations. The colors have the same meaning as in Fig. [3] 



where we assume that the change of mass provides a 
non- negligible contr i bution to the rate of p eriod change 
(jTurner et al.|[200l iNeilson fc Lester] 120081) . Since Po- 
laris has not been observed accreting material then the 
rate of change of stellar mass is due to mass loss only 
and is always < 0. This means that mass loss only acts 
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to increase the rate of period change. 

If we consider the largest value for the rate of pe- 
riod change during the blue loop evolution in Fig. [4] 
P w 1 s yr -1 , then to resolve the difference between the 
theoretical and observed rates of period change the rate 
of change of stellar mass must be < —7 x 10~ 6 Mq yr" 1 . 
We note this result ignores the possibility for feedback 
on the stellar radius due to mass loss, but the feedback 
is expected to decrease the mass-loss rate by less than an 
order of magnitude. 

5. DISCUSSION 

We hnd that the observed rate of period change, from 
all sources, is inconsistent with predictions from stel- 
lar evolution models unless one includes significant mass 
loss. Previous evidence for Cep heid mass loss is based on 
observations of infrared excess (INeilson et al.ll2009ll20lot 
iMarengo et al.ll2010t iBarmbv et al.ll201lD .~as well as 21- 
cm radio observations of a nebula surrounding 5 Cephei 
(jMatthews et al.ll2011[ ). These previous observations did 
not confirm mass loss because any circumstellar mate- 
rial could be relics from earlier stages of stellar evolu- 
tion or material coincidently existing along the observer's 
line-of-sight. Thus, the rate of period change provides a 
unique and strong argument for enhanced mass loss in 
Cepheids. 

The measured mass-loss rate for Polaris is simi- 
lar to that found from the 21cm radio ro observa- 
tions of 5 Cephei, as well as measurements from in- 
frared excesses, depend i ng on the a ssumed dust prop- 
erties. INeilson fc Lester! (|2008t 120091) predicted theoreti- 
cal mass-loss rates for Galactic Cepheids that are signif- 
icantly less than 10~ 6 Mq yr -1 , but stellar evolution 
models including the sam e pulsation-driven mass loss 
theory (jNeilson et aLl[~~0Tl~) found rates of the order 10 -7 
to 10 -6 Mq yr -1 . That work suggests that Cepheid mass 
loss may solve the Cepheid mass discrepancy assuming 
some moderate convective core overshooting. The signifi- 
cant amount of mass loss for Polaris implies that Cepheid 
pulsation is a wind driving mechanism. 

The measured rate of period change for Polaris is just 
one piece of evidence contributing to the understanding 
of Cephe id mass loss . Rec ent X-ray (jEvans et al.ll2010l l 
and UV (|Engle et al.l 120091 ) observations of Polaris indi- 
cate the presence of hot T ps 10 5 K plasma in the en- 
velope, which may origina te from shocks ge nerated by 
pulsation in the Cepheid (jFokin et al.l Il996| ) . Further 
spectral observations at near-UV wavelengths could pro- 
vide additional constraints on the mass loss mechanism 
in Cepheids. 

Polaris may not be unique as a Cepheid with sig- 
nificant mass loss. Large rates of period change are 
considered to be a property of first overton e pulsating 
Cepheids ()Szabadosl 11981 lEvans et al.ll2002j ). implying 
that many of these Cepheids, if not all, have significant 
mass-l oss rates during the blue loop stage of stellar evo- 
lution. I Turner et al.l (|2010[) found that four Cepheids are 
first overtone or double mode pulsators and argued they 
must be evolving along the Hcrtzsprung gap, based on 
the rates of period change. These are Polaris, DX Gem, 
BY Cas, and HDE 344787. However, we find that stellar 
evolution fits to the observed radius and effective tem- 
perature of Polaris suggest a rate of period change for 
the first crossing that is about four times larger than the 



observed rate; hence Polaris is most likely evolving along 
the blue loop. The observed rate of period change for 
Polaris is similar to that of BY Cas and DX Gem, while 
for HDE 344787, P = 12.96±2.41 s yr" 1 , suggesting that 
the rates for BY Cas and DX Gem may also be consis- 
tent with extreme mass loss for evolution along the blue 
loop. Observations constraining the radius, temperature, 
and/or luminosity along with rates of period change will 
test if all first overtone Cepheids have extreme mass loss. 

The results also suggest that more time-series observa- 
tions of Polaris in multiple wavelength regimes will help 
constrain the mass-loss history and mass-loss mechanism 
for this star and unravel a mystery of the North Star. 
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